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Abstract 

Optical spectroscopy in high magnetic fields B < 65 T is used to reveal the very different 
nature of earners in monolayer and bulk transition metal dichalcogenides. In monolayer WSe 2 , 
the exciton emission shifts linearly with the magnetic field and exhibits a splitting which orig¬ 
inates from the magnetic field induced valley splitting. The monolayer data can be described 
using a single particle picture with a Dirac-like Hamiltonian for massive Dirac fermions, with 
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an additional term to phenomenologically include the valley splitting. In contrast, in bulk 
WSe 2 where the inversion symmetry is restored, transmission measurements show a distinctly 
excitonic behavior with absorption to the Is and 2s states. Magnetic field induces a spin split¬ 
ting together with a small diamagnetic shift and cyclotron like behavior at high fields, which 
is best described within the hydrogen model. 

Keywords: Transition metal dichalcogenides, WSe2, monolayer, bulk, massive Dirac fermions, 
Fermi velocity, valley splitting 

Transition metal dichalcogenides (TX2, T=Mo,W,.. and X=S, Se,..) are quasi two dimensional 
layered materials with strong ionic-covalent bonding within a layer. A monolayer is composed 
of a single layer of transition metal embedded between two atomic layers of chalcogenides atoms 
in a trigonal prismatic structure. The weak van der Waals interlayer coupling nevertheless com¬ 
pletely modifies the band structure. Bulk crystals are indirect semiconductors with extremely weak 
phonon assisted photoluminescence (PL) emission. In monolayer TX2 the inversion symmetry is 
broken as the two sublattices are occupied by one T atom and two X2 atoms. The band gap is 
direct and located in two degenerate valleys (X±) at the corner of hexagonal Brillouin zone and 
the PL emission has a quantum yield which is four orders of magnitude larger than that of bulk 
crystals . 1 6 In addition, dielectric confinement, due to the very different dielectric environment 
outside of the monolayer, enhances the exciton binding energy beyond the usual factor of four 
limit for two-dimensional systems . 7 8 

Due to the d character of the orbitals, these materials exhibits a strong spin orbit coupling,® 
which when combined with the lack of inversion and time reversal symmetry, breaks the spin 
degeneracy and leads to a strong spin orbit induced Zeeman effect . 4 10 Hence, the spin and valley 
degrees of freedom are coupled which in turn leads to the valley contrasting selection rules for 
optical interband transitions with circularly polarized light : 11 12 o + couples to the K+ valley and 
a to K respectively (see inset in Fig [I]). Significant valley polarization was reported recently 
for monolayer M0S2 and WSe2 . 1315 Amongst the transition metal dichalcogenides, WSe2 has the 
most robust valley polarization, essentially due to the very large spin orbit splitting in the valence 
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band of this material . 13 

In transition metal dichalcogenides, the valley degree of freedom can be manipulated and po¬ 
tentially can be used to replace charge for the transmission of information . 11 13 15 19 For example, 
the anomalous Hall effect, whose sign depends on the valley index, was recently demonstrated 20 21 
and an external magnetic field can be used to tune the valley and spin polarization . 22 To date, the 
valley splitting has been investigated for WSe2 and MoSe2 in only moderately low magnetic fields 
B < 10 tJ 22 H !51 

In this letter we compare the optical response of monolayer and bulk WSe2 in high magnetic 
fields up to B — 65 T. We show that in monolayer WSe2 both the exciton and the trion exhibits 
a splitting which originates from the lifting of the valley degeneracy in a magnetic field. The 
magnetic field has little effect on the valley polarization of the neutral exciton, whereas it can 
significantly increase or decrease the valley polarization of the trion depending on the circular po¬ 
larization of the excitation light. The linear evolution of the energy of the exciton and trion features 
in magnetic field can be described using a Dirac-like Hamiltonian for massive Dirac fermions. In 
contrast, in bulk WSe2 where the inversion symmetry is restored, we see excitonic behavior with 
absorption to the Is and 2 s hydrogen like states. Both states show a spin (Zeeman) splitting accom¬ 
panied by a small diamagnetic shift at low magnetic fields before shifting as the cyclotron energy 
in the high field limit. 

For the measurements, single layer flakes of tungsten diselenide (WSe2) have been obtained 
by mechanical exfoliation of bulk 2 H-WSe 2 (the hexagonal 2 H-polytype of tungsten diselenide) 
single crystals grown using chemical vapor transport with Bromine as the transport agent. Samples 
obtained in this way are naturally n-type . 26 Two types of experiments have been performed; (i) 
micro-photoluminescence (/iPL) on direct gap monolayer WSe2 in steady state magnetic field 
up to 30 T and (ii) macro transmission of thin crystals of indirect gap bulk WSe2 in the pulsed 
magnetic field up to 65 T. Combining PL and transmission allows us to probe the A-exciton of an 
electron-hole pair at the K-points of the Brillouin zone in both monolayer and in bulk WSe2 for 
which PL is dominated by the indirect gap. 
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For /iPL measurements the WSe? monolayer flake was placed in a system composed of piezo¬ 
electric x — y — z translation stages and a microscope objective. The /iPL system was cooled in 
exchange gas to a temperature of T=4.2 K in a cryostat filled with liquid helium which was placed 
in a resistive magnet producing magnetic fields up to B — 30 T. The magnetic field was applied 
in the Faraday configuration and the sample was illuminated by a laser at 640 nm. The excitation 
power was kept low (of the order of a few hundred of nW) to avoid heating effects. Both the excit¬ 
ing and collected light were transmitted through a non polarizing cube beam splitter (50:50) placed 
on the optical axis of the objective. With this setup the circular polarization of the excitation and 
detection can be controlled independently. The diameter of the excitation beam on the sample was 
of the order of 1 /im. 

The magneto-transmission measurements were performed on thin bulk samples (also obtained 
by limited mechanical exfoliation) in pulsed fields < 65 T (~ 400 ms). A tungsten halogen lamp 
provides a broad spectrum in the visible and near infra-red range and the absorption is measured 
in the Faraday configuration with the c-axis of the sample parallel to magnetic field. Typical size 
of the spot was of the order of 200 /im and the polarization optics was introduced in situ. The 
sample was immersed in a helium bath cryostat pumped to a temperature of ~ 1.8 K to avoid 
noise due to the formation of helium gas bubbles. Note, that we are unable to perform macro 
transmission measurements on monolayer WSe 2 as only a tiny fraction of the transmitted light 
would pass through the micron size flake and our /iPL setup is not adapted for micro-transmission 
measurements. 

For both measurements the collected light was dispersed in a spectrometer equipped with a 
multichannel CCD camera cooled to liquid nitrogen temperatures. 

Fig[l]shows a comparison between the emission spectra measured on a single layer flake and the 
transmission through the thin bulk crystals of the WSe 2 all measured at zero magnetic field. In both 
spectra we observe a strong resonance around 1730 meV which corresponds well to the energy Is 
(n — 1) state of the A-exciton, previously reported for monolayer 13 27 29 and bulk WSe 2 - 27 30 The 
slightly higher energy of the A exciton emission observed in monolayer layer WSe 2 (compared to 
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Figure 1: Typical low temperature /iPL of monolayer and transmission spectra of bulk WSc 2 - 
The inset illustrates the optical valley selection rules in monolayer WSe 2 for o 1 and o circular 
polarization of the photons (c.b. labels the conduction band and v.b. labels the valence band). 
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bulk transmission) is also in agreement with previous studies. 27 28 We also see the 2s (n — 2) exci- 
ton level in transmission, previously reported in reflectance measurements. 30 The separation of the 
Is and 2s transitions, equal to 3i?*/4 in the 3D hydrogen model (E„ — R*/n 2 .n — 1,2,..) provides a 
rough estimate for the exciton binding energy in bulk WSe 2 of R* ~ 44 meV. The monolayer emis¬ 
sion spectrum has additional emission lines on the low energy side. The line around 30 meV lower 
than the exciton corresponds well to the negatively charged exciton (trion) emission, as our sam¬ 
ples are naturally n-doped. 26 Trion emission is routinely observed in mono layers dichalcogenides 
with an excess concentration of electrons and can be unambiguously identified by its characteristic 
power dependence. 31 32 Moreover, the observed binding energy of the trion corresponds well to the 
experimental value previously reported for WSe 2 - 13 29 3-1 Below the trion emission several peaks 
are observed which we label as “L” as they can be assigned to localized (bound exciton) states 
previously reported in WSe 2 and M 0 S 2. 15 29 33 

In the following we focus on the magnetic field dependence of the exciton transitions. Typical 
low temperature /iPL spectra of monolayer WSe 2 measured in dc magnetic fields up to 30 T are 
presented in Fig[2j a). The results of the transmission measurements through the thin layers of bulk 
WSe 2 measured in pulsed magnetic fields up to 65 T are presented in Fig [2f b). For both data sets 
the solid and dotted lines correspond to o + and o detection respectively. The photoluminescence 
was excited using o + polarization while for the transmission unpolarized white light was used for 
the excitation. 

In both bulk and monolayer crystals the polarization resolved measurements reveal a splitting 
of the exciton transition, nevertheless with very different physical origins due to the markedly dif¬ 
ferent optical selection rules in bulk and monolayer WSe 2 - Additionally, the monolayer emission 
shows that, already at zero magnetic field, both the neutral and charged exciton populations are 
partially valley polarized as the emission intensity is not the same for o and <7 detection. This 
is not surprising, as our circular excitation was relatively close to the exciton resonance, as also 
previously reported for WSc 2 -' 13 

In monolayer WSe 2 the circular polarization of the excitation allows us to selectively pump 
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Figure 2: Typical low temperature spectra at various magnetic fields of (a) /iPL from monolayer 
WSe 2 and (b) transmission through the thin WSe 2 bulk crystal, using a + and a~ detection. For 
the PL shown here <7 + polarized light was used for excitation while the white light for the trans¬ 
mission measurements was unpolarized. The dashed lines are drawn as a guide for eye showing 
the evolution of exciton features with magnetic field. The spectra have been shifted vertically for 
clarity. 
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Figure 3: (a),(b) Typical /IPL spectra excited o ! and o polarized laser respectively. The red/blue 
lines denote the detection in o 1 and a respectively. (c),(d) degree of circular polarization of the 
neutral and charged exciton respectively. The open symbols denote o + excitation as closed one 
correspond to the o polarization of the excitation. 


the K + or K valley. The circular polarization of the emission then depends on the inter valley 
scattering rate and the radiative recombination lifetime. In the absence of inter valley scattering 
the emission would be 100% polarized. The degree of circular polarization P for the exciton and 
trion emission can be defined as 


P 


Ax+ At 
V+At-’ 


( 1 ) 


where I 0 ± is the intensity of the photoluminescence for a given circular (detection) polarization 
obtained by fitting the spectra with a Gaussian function for each value of magnetic field. 

We have analyzed data for all four possible combinations of the circular polarization of the 
excitation and detection. Typical polarization resolved /iPL spectra measured at 30 T for cr + and 
(7 are presented in Fig[3ja) and (b) respectively. Note that the direction of the magnetic field was 
the same for all measurements. The required polarization combination was selected by rotating 
the excitation and detection polarization optics located outside of the cryostat. As in zero field, the 
30 T raw data clearly shows that both neutral and charge exciton exhibits a significant degree of 
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circular polarization. 

The full magnetic field dependence of calculated degree of circular polarization of the neutral 
and charged exciton calculated using Equation 1 is presented in Fig |3fc),(d) respectively. Both 
neutral and charged exciton shows a valley polarization in zero magnetic field. For the exciton 
P = 25% and for the charged exciton P = 40%. Although we are not exciting exactly at the res¬ 
onance it has been shown that WSe 2 , which has the largest spin orbit induced spin splitting in 
both valence and the conduction bands of any member of the TX 2 family, has a valley polariza¬ 
tion which is the most robust of all the transition metal dichalcogcnidcs. 11 For the exciton the 
degree of the circular polarization remains almost unchanged with the magnetic field suggesting 
that magnetic field has little influence on the inter valley exciton scattering rate. On the other hand 
the degree of polarization is trion is strongly influenced with the polarization of the (J 1 excited 
emission (pumping K+ valley) increasing strongly, accompanied by an equally strong reduction in 
the degree of polarization for o excitation (pumping K- valley). This suggests, that although the 
magnetic field does not modify the inter valley exciton scattering rate, it can nevertheless modify 
the probability for trion formation, most likely via its influence on the excess electron population 
of the spin orbit split spin levels in the conduction band for each valley. 

The evolution of the exciton and trion polarization reported here is somewhat different from 
previous investigations of WSc 2^ 2 and MoSc 2^ 3 which might be due to the low magnetic fields 
employed (B < 7 T). For example, in our data, the true behavior of the trion polarization is only 
revealed for magnetic fields B > 5 T. Additionally, the polarization depends sensitively upon the 
exact excitation used, together with sample dependent parameters such as the valley scattering rate. 

The energy of the polarization resolved exciton and trion emission for monolayer WSe 2 is 
shown in figure [4j a). Both the exciton and trion energy shifts linearly to higher energy with in¬ 
creasing magnetic field showing a pronounced splitting which also evolves linearly with magnetic 
field. In the transition metal dichalcogenides the carriers behave as massive Dirac fermions which 
can be described by a Dirac like Hamiltonian. 11 34 In this simple picture the magnetic field quan- 
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tizes the energy of the carriers in the conduction and valence bands into Landau levels with energy, 


E^ k = A \/A 2 + ne 2 ± ^g vJ u 5 fl, (2) 

where A = ± 1 designates the conduction/valence band, A is half of the band gap and £ = \flhvp /(« 
is the characteristic magnetic energy. Here n is the Landau level index, vy is the Fermi velocity and 
£g — (h/eB) is the magnetic length. Note in this single particle picture 2A is obtained from the 
observed zero magnetic field transition energy to which the exciton binding energy in monolayer 
WSe 2 of Eb — 370 meV has to be added. 35 



Figure 4: (a) Magnetic field dependence of the energy of the polarization resolved exciton emis¬ 
sion in monolayer WSe 2 measured at T — 4.2 K. The solid lines are fits to the data using the 
Hamiltonian for massive Dirac fermions, (b) Magnetic field dependence of transmission in bulk 
WSe 2 measured at T = 1.8 K. The solid lines are fits to the data using Equation 3 to model the 
diamagnetic (quadratic) shift at low magnetic fields followed by a linear Landau like behavior at 
high magnetic fields. 


As in graphene, in the transition metal dichalcogenides the selections rules for dipole allowed 
optical transitions are An — ±1, i.e. the Landau level index changes by one. In this case spin is 
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conserved and any spin splitting of the Landau levels will not give rise to a splitting of the optical 
transition unless the electron and hole g-factors are significantly different. Thus, the observed 
splitting of the transitions is included simply as a phenomenological valley splitting , /!/;/:> 
where g v is the effective valley g-factor. The n > 0 Landau levels always occur in pairs with one 
in each valley in both the conduction and valence bands. The n = 0 level is special and there is a 
single n = 0 Landau level per valley. For both spin levels, the n = 0 Landau level is fixed at the 
top of the valence band for the +K valley and at the bottom of the conduction band for the —K 
valley. 34 Thus, in a magnetic field resonant c ± excitation selects 0 —H or 1 —> 0 optical transitions 
and therefore selecting the valley as in the zero magnetic field case (see schematic in Figure [5]). 


K+ K- 



Figure 5: Schematic showing evolution of the n — 0,1 Landau levels of massive Dirac fermions in a 
magnetic field together with the dipole allowed 0 —> 1 and 1 —> 0 transitions which select the valley 
involved (solid lines are Landau levels, solid arrows indicate the transitions). Spin is conserved so 
that the Zeeman energy (not shown here) does not split the transition. Any asymmetry of the 
electron/hole Dirac cones, which modifies the energy of the n — 1 Landau levels, would split the 
transition as indicated by the dashed lines/arrows. 


The solid lines in figure [4j a) are the transitions energies calculated from the energies of the 0 —> 
1 and 1 —)• 0 transitions from Equation 2 with a Fermi velocity vp = (0.51 ±0.02) x 10 6 m/s and an 
effective valley g-factor g v — 4 ±0.5. We have included excitonic effects simply by subtracting the 
exciton (trion) binding energy from the calculated transition energy. The parameter set provides a 
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reasonable fit to both the exciton and trion emission data. The transition energies evolve linearly 
with magnetic field which can be understood as follows: The band gap of WSe 2 is large so that 
for the magnetic fields of interest here A 2 £ 2 so that the orbital contribution to the Landau level 
energy is well approximated by a Taylor expansion 

E n — ±A(1 +/7£ 2 /A 2 )2 ~ ±(A-\- ntieBx’]? / A). 

Thus, the carriers behave as massive Dirac fermions with an effective mass m* — A/Vp ~ 0.1 m e 
where m e is the free electron mass. Note, that this result can also be obtained by comparison with 
the Klein-Gordon equation E 2 = m 2 c 4 + p 2 c 2 replacing the speed of light with the Fermi velocity. 
For this reason A = me 2 is often referred to as the “mass gap”. 

The observed splitting is asymmetric due to the orbital (Landau level) contribution which is 
always positive. The low energy peak almost does not move in magnetic field as the orbital and 
valley contributions are of similar size but of opposite sign. Li et al. 25 observe a similar valley 
splitting in MoSe 2 of 0.24meV/T but with a symmetric splitting in fields up to 10T. This may 
suggest that the orbital contribution is quenched in their measurements, either due to the higher 
temperature used (T = 10 K) or increased disorder in their sample (ct) c T < 1). 

In the literature the splitting of the transition is often described, using slightly different lan¬ 
guage, in terms of a two band tight binding model in which the atomic J-orbitals which form the 
valence band have an in plane angular momentum mh with in — ±2 according to the valley index. 
In magnetic field, the valley splitting is then 4/igfi which agrees within error with the splitting ob¬ 
served here (we estimate from the fit that g v — 4.0 ± 0.5). From low field measurements (B < 107) 
Srivastava and co-workers find a similar value for the valley splitting of WSe 2 (g v = 4.3)J 24 while a 
factor of two lower value (~ 2 BbB) has been the published recently by Aivazian and co-workers. 22 
In the tight binding model, the valley magnetic moment due to the self rotation of the wave packet 
gives no contribution to the splitting provided the electron and hole masses are identical. Aivazian 
et al. 22 invoked different electron and hole masses to correctly predict the observed splitting in 
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their samples. Similarly, in our Hamiltonian any asymmetry of the electron and hole Dirac cones 
would lead to a splitting of the 0 —>■ 1 and 1 —> 0 transitions (see Figure [5]). However, as our data 
corresponds within error to the expected splitting of 4/!«/:>. this suggest that electron-hole asym¬ 
metry in WSe 2 is small and makes no significant contribution to the observed valley splitting. 

For bulk WSe 2 , the dipole selection rules are An = 0 so that the angular momentum of the 
photon is used to flip the spin. The energy of the polarization resolved exciton features in trans¬ 
mission for bulk WSe 2 shown in figure [4} b) clearly shows signs of excitonic effects with quadratic 
diamagnetic shift at low magnetic fields which becomes more linear at high fields. The evolution 
of the data in magnetic field can be well described using 

, ^ ft(Oo h / o o 1 

=-^+2V ft) o + ® c ± 2 Ss ^ bB ' (3) 


where Hcoq is an energy which controls the diamagnetic shift, h(D c = eB/m* is the cyclotron energy 
and m* is the reduced exciton mass. The observed spin splitting is included via an effective g-factor 

8s- 

Neglecting spin for the moment, at high magnetic field ft ) 2 Q)q and the energy is well ap¬ 
proximated by hco c /2. At low fields when (Dq A> ft) 2 to first order a Taylor expansion gives the 
diamagnetic shift 

he 2 B 2 ? ? ? 

A E ~ -- T = e~anB~. 

4 (Dom* 2 u 

Here we have defined a length ao — ^ /h/4m* 2 (Oo which can be identified with the effective Bohr 
radius in the hydrogen model. 

The solid lines in figure |4jb) are the transition energies calculated using Equation 3. For the Is 
absorption, we use an effective mass m* — 0.1m e as in monolayer WSe 2 - The spin splitting can be 
reproduced using an effective g-factor g s = 2.3 and the diamagnetic shift reproduced using Hcoq — 
8 meV. This gives an estimation of the diamagnetic shift at low magnetic fields of A E « 0.9/iV/T 2 . 
The effective Bohr radius is ao ~ 1.8 nm from which we can obtain a rough estimate of the exciton 
binding energy in the hydrogen model e 2 /%TZ£ r £oao « 56 meV using a dielectric constant £ r ~ 7 for 
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bulk 2H-WSe2. This value agrees well with the exciton binding energy of ~ 55 meV determined 
directly from the n = 1,2,3 hydrogen series of levels in reflectance measurements on bulk 2H- 
WSe 2 crystals. 30 A similar binding energy (~ 56 meV) has been reported 36 in high magnetic field 
measurements up to 150 T on bulk M 0 S 2 , however with a much smaller diamagnetic shift ~ 
0.2/ieV/T 2 . The approximately four fold smaller diamagnetic shift is however, fully consistent 
with our data on WSe 2 when the smaller Bohr radius (1.28 nm) and smaller reduced exciton mass 
(0.4m e ) in M 0 S 2 are taken into account. 36 

For completeness, we have also fitted the 2s data using the same effective mass in* — 0.1m e , 
a slightly larger effective g-factor is required g s — 2.8, with a smaller value of ha>o = 2 meV. The 
estimated diamagnetic shift A E « 3.4/iV/T 2 is larger due to the more delocalized nature of the 2s 
wave function which extends over « 3.7 nm. 

In conclusion, we have investigated monolayer and bulk WSe 2 in high magnetic fields up to 
B = 65 T using optical spectroscopy. In monolayer WSe 2 , the exciton emission exhibits a splitting 
which originates from lifting of the valley degeneracy by the magnetic field. The linear evolution of 
the energy of the exciton features in magnetic field can be described using a Dirac-like Hamiltonian 
for massive Dirac fermions with a Fermi velocity of 0.51 x 10 6 m/s and an effective valley g-factor 
g v — 4 for an assumed effective spin 5=1/2 system. The measured Fermi velocity can be used 
to estimate the effective hopping integral t of the tight binding Hamiltonian using at — hvp where 
a = 3.31. The measured value of hvp — 3.36eV suggests t — 1.02 eV which is roughly 15% less 
than the effective hopping integral predicted from first principle band structure calculations. 11 34 
In contrast, in bulk WSe 2 where the inversion symmetry is restored, transmission measurements 
show that the exciton exhibits a spin (Zeeman) splitting and the exciton Is and 2s features show a 
small diamagnetic shift which allows us to determine the exciton binding energy of around 56 meV 
within the three dimensional hydrogen model. In two dimensions, the binding energy is enhanced 
by at most a factor of four, so the predicted 2D exciton binding energy (4 x 56 ~ 224 meV) is 
significantly less than the established value of 370 meV in monolayer WSe 2. 30 This suggests that 
dielectric screening (image charge), due to the very different dielectric environment outside of 
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the material,— significantly enhances the exciton binding energy in monolayer transition metal 
dichalcogenides. 
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